Resonance parameters of the rho-meson from asymmetrical lattices 
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We present a lattice QCD calculation of the parameters of the p meson decay. The study is carried 
out on spatially asymmetric boxes using nHYP-smeared clover fermions with two mass-degenerate 
quark flavors. Our calculations are carried out at a pion mass m v = 304(2) MeV on the set of 
lattices V = 24 2 x t?24 x 48 with r\ = 1.0, 1.25, and 2.0 with lattice spacing a = 0.1255(7) fm. The 
resonance mass m p = 827(3) (5) MeV and coupling constant g P i„r — 6.67(42) are calculate using 
the P-wave scattering phase shifts. We construct a 2x2 correlation matrix to extract the energy 
of the scattering states and compute the phase shifts using the finite volume formula. By varying 
the degree of asymmetry, we are able to compute a set of phase shifts that are evenly distributed 
throughout the spectral region where the p decays. 

PACS numbers: 11.15.Ha,12.38.Gc 



I. INTRODUCTION 

The study of multi-particle systems is an important 
step for lattice QCD and our understanding of the strong 
force. Experimentally, the interaction of elementary par- 
ticles is studied using scattering methods. Due to bet- 
ter algorithms, more efficient codes, and an increase in 
computational resources, it is finally possible to use lat- 
tice QCD to compute scattering observables. Of par- 
ticular interest are the scattering channels that exhibit 
resonances. In this work, we focus on the light vector 
meson p(770) seen as a resonance in the elastic scatter- 
ing of two-pions in the I G (J PC ) = 1 + (1~~) channel. 

Currently, the method of choice to carry out lattice 
studies of scattering phase shifts is the Lusher method [1- 
4]. Lusher showed that scattering phase shifts can be 
computed from the spectrum of two-particle states on a 
torus. This method circumvents the "Maiani- Testa no- 
go theorem" for Euclidean field theories [5]. To mea- 
sure scattering phase shifts then requires the computa- 
tion of the energies of the two-particle states. While it 
is known how to compute the energies of the scattering 
states, the task can be quite challenging and is computa- 
tionally costly. In the case of the p decay, as the physical 
pion mass is approached the scattering states get closer 
together and we have to carefully disentangle the relevant 
two-pion states. Additionally, when the four-pion states 
become dynamically relevant, the formula is no longer 
valid. 

Lusher's formula was developed for lattices with cu- 
bic symmetry. Two notable extensions were worked out: 
(1) systems with non-zero total momentum where the cu- 
bic symmetry is lost due to the rclativistic boost to the 
center-of-mass frame [6] and (2) systems with zero to- 
tal momentum on asymmetrical lattices [7, 8]. The first 
method can be used to investigate the p decay on lattices 
where the zero-momentum p cannot decay due to kine- 
matical constraints — this was used recently to study the 
p resonance [9-11]. Using this method different phase 
shifts can be extracted by varying the total momentum 
of the system. While it has been successful for studying 



the p resonance, the ability to tune the momentum of the 
two-pion system is limited by the finite step size of the 
momenta. 

In this work we chose to use the second extension. This 
method offers a finer control over the two-pion momen- 
tum. With a finer control over the momentum, it is pos- 
sible to map out narrower resonances, and as will be dis- 
cussed later, the overall increase in computational cost 
is not that significant. Another reason for using asym- 
metrical lattices is that the projection onto the relevant 
irreducible representation can be performed cleanly. This 
projection is always performed in the center-of-mass, 
which coincides with the lattice frame when the interpo- 
lators have zero momentum. For system with non-zero 
momentum, the boost mixes representations and we have 
to rely on dynamics to disentangle the states of interest. 

In this work we employ nHYP-smeared clover fermions 
with two mass-degenerate quarks. The simulation is car- 
ried out with a pion mass = 304(2) MeV at a lattice 
spacing a = 0.1255(7) fm. Wc construct a 2 x 2 cor- 
relation matrix and extract the ground and first excited 
state energies on dynamically generated ensembles of 300 
gauge configurations for the three lattice volumes 



V = 24 2 x r?24 x 48 , rj = 1.0, 1.25, 2.0 



(1) 



We extract six scattering phase shifts well distributed 
throughout the resonance region and compute the reso- 
nance mass and coupling constant. The methods used 
here were first tested on a set of quenched configurations 
for similar quark masses and lattice sizes [12]. 

The paper is organized as follows. In Section II we 
describe the methods used and collect the relevant for- 
mulas. In Section III wc present our results for the two- 
pion spectrum, compute the phase shifts, extract the res- 
onance parameters, and compare our results with other 
studies. In Section IV we present our conclusions and 
discuss future plans. 
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II. METHODS 



Phase shift formulas 



TABLE I: Resolution of the 2J + 1 spherical harmonics into 
the irreducible representations of Oh and D^. 



To compute the elastic scattering phase shifts, we use 
the well known Liisher's formula [1-4]. In particular, we 
make use of its extension to asymmetrical lattices [7, 8]. 
In the following we list the relevant formula needed to 
determine the scattering phase shifts 5i(k) for two pions 
with back-to-back momentum in the angular momentum 
I = 1 channel. 

In this work we consider lattices with one spatial di- 
rection elongated. Since we work on a spatial torus, the 
symmetry group for zero-momentum states is reduced 
from SO(3) to a discrete subgroup. In this case the rel- 
evant symmetry group is D±h- For concreteness we take 
the z-direction to be elongated. The decomposition of the 
irreducible representations of SO(3) into D±h are listed 
Table I. From the table we see that Z = 1 is the low- 
est angular momentum mode which couples to the 
and E~ representations; therefore, either representation 
is suitable. In the following we present the formula for 
these two channels. 

Following [8], we introduce the generalized zeta func- 
tions 



O h 



Zi m (s,q 2 ;r] 2 ,m) = ^2 



3Wn) 
(n 2 - q 2 ) s 



where yi m {i\) are the harmonic polynomials 
y h n(n) = h l Yi m (n & ) 

and 

n = (ni,n 2 /?72, "3/773), neZ. 



(2) 



(3) 



(4) 



The series in Eq. (2) is convergent for Re 2s > I + 3 and 
can be analytically continued to the half plane Re 2s > 
1/2. To compute the scattering phase shifts, we need to 
evaluate Zi m (s = 1, q 2 ; 772, 773). The details are given in 
Appendix A. The parameter q is related to the invariant 
energy W of the two-pion system through the relation 



kL 
2^ 



(5) 

where fc, the pion "momentum" , is determined by 

W = + k 2 . (6) 

The 3D torus geometry and spatial volume is V = L x 
772L x r] 3 L. For notational simplicity, we also introduce 
the quantity 



Wz m (l,q 2 ; 772,773) 



Z; m (l,q 2 ; 772 , 77 3 ) 
7T 3/2 772773g' +1 



(7) 



Since our lattices are elongated only in one direction, we 
set 772 = 1 and 773 = 77 > 1. The phase shift formula for 
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the two representations are then given by 



A 2 : cot St (k) = Woo + -j=W 2 o 

V5 

E-: cot 6i(k)= Woo -4? W 20 

V5 



(8) 
(9) 



where contributions from angular momenta / > 3 are 
assumed to be negligible. 

We note that our formulation is slightly different 
than [8], which defines the lattice volume as V = r/iL x 
1J2L x L. The phase shift formulas in either case are equiv- 
alent. Additionally, our formula for the E ~ channel dif- 
fers by the exclusion of the term ±^/3/10(W22 + W2-2), 
which vanishes due to symmetry under rotations around 
the elongated direction. 

Lastly, we consider the case of a cubic lattice. The 
symmetry group is Oh, and the relevant representation 
is Fi (sec Table I). In this case W20 vanishes, Eq. (8) and 
Eq. (9) become equivalent, and the phase shift formula 
is given by setting the right hand side equal to Woo- 



B. Selecting lattices 

In order to extract the resonance parameters for p, we 
need to compute the scattering phase shifts in the spec- 
tral region where the resonance appears. To minimize 
statistical errors, it is preferable to compute the phase 
shifts using the lowest lying states in the channel. To 
achieve this, we have to select appropriate lattice vol- 
umes and, for practical purposes, make a selection which 
is computationally economic. This choice is difficult since 
the resonance mass and coupling constant are generally 
not known a priori for unphysical pion masses. However, 
using recent lattice results [9, 11, 13] and predictions from 
unitarized chiral perturbation theory [14], a reasonably 
good guess can be made. In the following we describe 
the method we used to select lattices. 

To estimate the two-pion energies, we use the effective 
range formula [15] 



cot 61 (WO = 



6tt W(m 2 p -W 2 ) 



(10) 
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which is known to parameterized the p resonance well 
at the physical pion mass. Using the equation above 
together with Eqs. (8) and (9), we can compute the ex- 
pected spectrum as a function of 77 provided we can es- 
timate the values of m p and <7 p7 r7r- From recent lattice 
studies we know that the coupling constant has little 
pion mass dependence and is close to the physical value 
g p7r ^ m 6 for pion masses < 350 MeV [9, 11, 13]. As 
a result, we start with the assumption that <7 pw7r = 6. 
To estimate the value of m p , we use the result from uni- 
tarized chiral perturbation theory which shows reason- 
able agreement with the lattice results [14]. Note that 
this picture can be refined as we collect data from lattice 
simulations. When choosing subsequent lattices, we use 
estimates for m p and g p7r7T based on the available lattice 
data. 

In Fig. 1 we plot the expected spectrum for a pion mass 
of 300 MeV. Note that the two-dimensional E~ represen- 
tation corresponds to the two states with non-zero back- 
to-back momentum in the transversal directions, whereas 
the A2 is the one-dimensional representation correspond- 
ing to the state with the momentum along the longitudi- 
nal/elongated direction. For the E~ representation the 
spectrum depends very little on 77. The reason is that the 
relative momentum of the two-pion state depends on the 
size of the transversal direction, which does not change. 
Thus, computing the energy of the states in the E~ rep- 
resentation will produce phase shifts very similar to the 
ones on a cubic lattice. The representation shows a 
strong dependence on the value of 77 and offers good con- 
trol over the value of the ground and first excited state 
energies. To generate additional phase shifts, one then 
needs to work in the channel and select appropriate 
values of r\. We see that by choosing 77 = 1, 1.25, and 
2.0, we should get a good scan of the resonance region. 

Lastly, we make a few remarks about our choice to use 
lattices with a single direction elongated. It may seem 
that a better choice for these types of studies are or- 
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FIG. 1: Estimated spectrum for — 300 MeV, <? p7r ,r = 6, 
and m p — 830 MeV. The lattice size is set to L — 3 fm. The 
two curves indicate the ground, and first excited states. The 
shaded region is the resonance region, i.e., m p ± T p . 



thorhombic lattices, i.e., lattices in which all three side 
have different lengths. In this case six different phase 
shifts can be computed from one ensemble using only the 
ground and first excited state energies. Choosing appro- 
priate values for L,r]2, and 773, we could get all of these 
phase shifts in the resonance region. Since the computa- 
tional effort scales with V 5 / 4 [16], in this case the compu- 
tational cost is (rj2T]3) 5 / 4 x r where r denotes the cost for 
the cubic lattice. If instead we choose a lattice with only 
one direction elongated and generate the two ensembles 
with 77 = 7/2 and 77 = 773, the cost is (r^J 4 + ??3^ 4 ) X r. The 
overall computational cost of the two choices is similar 
when the elongation parameters 772,3 are small. Addition- 
ally, using two different ensembles reduces the correlation 
between phase shifts and allows for studies of finite vol- 
ume effects of other observables. In light of these consid- 
erations, we choose to generate lattices with only a single 
direction elongated. 

C. Variational analysis 

To compute the low lying two-pion energies with the 
quantum numbers I G (J PC ) = 1 + (1 ) of the p 1 we em- 
ploy the variational method proposed by Liischer and 
Wolff [17]. We use a two-dimensional variational basis 
and construct the correlation matrix 

C(% = <0<(*K>J(O)) (11) 

to extract the ground and first excited states. We com- 
pute the eigenvalues of 

C(t )- 1/2 C(t)C(7io)- 1/2 V' (n) (t,*o) = A (n) (Mo)l& (w) (*,*>) 

(12) 

for each time slice t. The two-pion energies are then 
determined from the long-time behavior of the eigenval- 
ues [18] 

A (,l) (i,io)oce- B " t (l + O(e- Af; " t )), 71= 1,2 (13) 

where AE n = E3 — E n . For the variational basis we use 
the interpolators 

01 = ^{7r+(p>-(-pJ-7r-(pj7r+(-pJ}(,14) 

02 = ^-^={u(x) l3 u(x)-d(x) l3 d(x)} (15) 

where p z = is the lowest lattice momentum in the 
elongated direction and 

7T-(p) = ^«(x) 75 d(x)e i P x , (16) 

X 

7r+(p) = ^J(x) 75U (x)e*P x . (17) 

x 

Lastly, we note that as defined above Oi and O2 are 
hermitian and anti-hermitian, respectively 
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D. Evaluation of the correlation matrix 

To evaluate the correlation matrix, one needs to com- 
pute the following diagrams: 

Cu(t) 

□ □ X X U 00 



c 12 (t) 



C 21 (t) 



c 22 (t) 



VV A A 

Above, we used a compressed notation; more explicitly, 
we have 

(-p,t/) (p,t f ) 

75 r-* 1 75 

I | = Tt[e^iiM-\t i ,U)e ip '1sM- 1 (t i ,tf)x 

(-l 5 M) (pi) e*nM- 1 (t f ,t f )e- i ''T B M-\t f ,t i )], («) 

where the square matrices 

(e lp ) x ,a, a ; V ,b,(3 = e ipx S(x - y)8 a b5 a p , (19) 
M~ 1 (tf,t i ) x , a , a]y , b ^ = M~ 1 (x,tf,a,a;y,t i ,b,l3), 

have 12 x N s rows with N s the number of points in a 
time slice. Above, M is the quark matrix and thus M _1 
is the regular quark propagator. All four point diagrams 
arc defined in a similar manner. The diagrams with three 
points can be written using the following template: 



(o,t f ) 

73 



Tr [e- ip ^ 5 M-\u,ti) x 



75 75 

i-P,n) (p,u) e ip 75M- 1 (^,t / )7 3 M- 1 (t/,i i )]. (20) 

Using the 75-hermiticity of the quark propagator, i.e., 
M _1 (t, t'Y = 7 5 M _1 (t',<)75, we can show that the four- 
point diagrams are purely real and the three-point dia- 
grams arc purely imaginary, configuration by configura- 
tion. Moreover, using parity we can show that 



□>=<□>„• 



V -- V 



(21) 



Time reversal symmetry allows us to connect the quark 
diagrams required for C\ 2 with the ones for C 2 \. In par- 
ticular, one can show that 



(22) 



Above, we use to denote the average with respect to 
the gauge configurations. All together, the four compo- 
nents of the correlation matrix can be constructed from 



- 2 



+ 



Gu(t) = {2^\ 

c 12 (t) = -c 2 i(t) = , 
c 22 (t) = 



00 



(23) 



Note that C\ 2 and C 2 \ are related. To save time, we 
compute only C 2 \ and use it to determine G\ 2 . C 22 can 
be evaluated using standard methods but the terms for 
Cxi and C21 require the all-to-all propagator. As a result, 
they must be computed stochastically. Our calculation 
follows the steps described in a study by the CP-PACS 
collaboration [11]. We introduce a random Z(A) noise 
time-slice vector £ (with color and spinor components), 
satisfying the condition 



<«*> 



S 



1 



(24) 



where 1 is the identity matrix in 12 x N s dimensions and 
denotes the average over the noise. Using these noise 
vectors, we define the time-slice vectors 

u(t\p,t 2 ,0 = M- l {t,t 2 y p Z (25) 
«(t|pi,ti|p2,*2,0 = Af- 1 (t,i 1 )e <l ' I 7feu(ti|p,t a ,0- 

With the above definitions unbiased estimators for the 
quark diagrams for C\\(t) and C 2 ±(t) can be constructed 
as 



□ 



(v(t f \~p : U\p,U,0^e- tp v{t f \p,U\0,U,0)^ u , 

(tt(t/|0,t i ,0 + e~ <p tt(t/|p,ti,0>£ x 
{u(t f \0,U,7 i ye lp u(t f \~p 1 U,7 1 )) r 



(26) 



00 = (<"(*/|0,ti,6 + e-^(t/|-P,*i,e)> s x 



A 
A 



(u(t f \0,t i ,r 1 ye lp u(t f \p,t i ,r ] )) ri 
(u(t f \0, U, £) + 7s73«(*/l - P, U\p, U, 0)^ v 
(u(t f \0, U, t 7573«(i/|p, U\ - P, k, 0) LU 



Lastly, we make a few comments on the variance of 
the stochastic estimators. As is common practice, we 
chose to use spin-color dilution. To check its effective- 
ness, we compared the variance of Cn(t) and C 2 \(t) with 
and without dilution in the temporal range where the fit- 
ting is performed. Using spin-color dilution, we found an 
order of magnitude decrease in the variance for the same 
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FIG. 2: Left panel: ground state energies for the three different elongations extracted using different methods. Right panel: 
ground and first excited state energies compared to the effective range formula using the parameters extracted from our fits. 
The dashed lines represent the expected two-pion energies in the absence of interaction, i.e., W„ = 2\Jm\ + {2-Kn/rjL) 2 . 



numerical effort. Another decision which has to be made 
is the size of the noise ensembles. Using a simple model, 
we can show that the overall variance of the estimator is 

^#(' + ^) 

where ajj and tr| are the contributions to the overall vari- 
ance from the gauge and noise fluctuations, N is the 
number of gauge configurations, and M the number of 
noises per configuration. In the temporal region where 
the fits were performed, we found that the ratio ct|/<7^- is 
small and that the variance is dominated by gauge fluc- 
tuations. As a result, the variance is not significantly 
decreased even after a single noise. This is similar to the 
results found by [11] where no substantial decrease in the 
variance was found after two noises. In this work, we also 
chose to use two noises. 

The stochastic estimator requires thousands of inver- 
sions for each configuration making the numerical cost 
of the analysis comparable with the cost of generating 
the dynamical configurations. To carry out these calcu- 
lations, we used our own GPU BiCGstab inverters [19]. 

III. RESULTS 



with temporal extent N t = 48. The lattice spacing for 
all three ensembles was determined to be a = 0.1255(7) 
using the Sommer scale (see Appendix B), and the com- 
puted value for the pion mass was = 304(2) MeV. 

We only computed the energies in the channel since 
the energies in the E~ channel are not significantly dif- 
ferent than those obtained in the F~ channel on the cu- 
bic lattice. The numerical values are listed in Table II. 
In the left panel of Fig. 2, we show the ground state 
energies obtained using solely the 7T7T and p interpola- 
tors, as well as the energy obtained using the variational 
method. The results show consistent behavior as the lat- 
tice size is increased with the energy and, as expected, 
decreases as the back-to-back momentum of the pions be- 
comes smaller. Note that the individual correlators, ex- 
tracted separately from Cu(t) and C*22(i), have larger er- 
ror bars. For our fits we use the energies extracted using 
the variational method. In the right panel of Fig. 2, we 
show our results for the ground and first excited states, 
as well as the spectral behavior as predicted by the ef- 
fective range formula Eq. (10) using our fit results. We 
see that the effective range formula describes the lattice 
data well and that the lattice data resolves clearly the 
effect of the pions' interaction. 



A. Two-pion spectrum 



To compute the two-pion spectrum, we generated three 
dynamical ensembles of 300 gauge configurations. We 
employed the Lusher- Weiss action with f3 = 7.1, re = 
0.1282, and uq = 0.868. The simulation was carried 
out using nHYP-smeared clover fermions with two mass- 
degenerated quarks and the standard smearing parame- 
ters ai = 0.75 ,a2 = 0.6 , and 013 = 0.3 [20]. The spatial 
volume of the the three lattices was 

24 x 24 x 7^24, 77 = 1.0 , 1.25, 2.0 , (28) 



TABLE II: Two pion energies in then A 2 sector. 



V 


1.0 


1.25 


2.0 


aEo 


0.516(10) 


0.511(4) 


0.442(4) 


aE 1 


0.660(6) 


0.606(10) 


0.549(12) 


Si(E ) 


1.71(11) 


0.822(65) 


0.242(31) 


Si(Ei) 


2.99(10) 


2.64(13) 


2.04(18) 


am-,, 


0.1925(7) 


0.1944(6) 


0.1946(8) 
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FIG. 3: Left panel: scattering phase shifts computed in this study, effective range formula fit, and the results from unitarized 
XPT [14]. The dashed lines describe the energy dependence of the phase shift according to Eq. (8). Right panel: cross section 
as a function of the center-of-mass energy. The dashed line indicates the p mass, the narrow dark band its error, and the 
shadowed box indicates the resonance region m p ± T. 



B. Phase shifts and resonance parameters 

To compute the resonance parameters, we evaluate the 
scattering formula Eq. (8) with the energies obtained for 
the ground and first excited states. The resonance mass 
and coupling constant arc then computed by fitting the 
scattering phase shifts with the effective range formula 
given in Eq. (10). We use a correlated x 2_n t since the en- 
ergies for the ground and first excited states are extracted 
from the same ensemble and construct the covariance ma- 
trix with two-by-two blocks, one per ensemble. 

Note that the standard x 2 method needs to be modi- 
fied for our fit: the statistical errors affect both x and y 
coordinates in our plot since the phase shifts are related 
to the energy. Referring to the left panel of Fig. 3, as the 
energy is varied the phase shifts move along the dashed 
lines indicated in the figure. The solid lines indicate the 
error bars obtained by varying the energy in the W ± aw 
interval. The \ 2 function is defined to be 



1 



X 



A T C _1 A 



where A is the difference vector 



A, = Wi - W° 



(29) 



(30) 



Wi indicates the computed value, and W® the value at 
the intersection of the dashed line and the fit curve. 

The decay width is determined using the fit results and 
the relation 



9%* ( m2 P /^ ~ ™l) 
6tt ml 



3/2 



We find 
g pir7r = 6.67(42) . 



m p = 827(3) (5) MeV, 



T = 76.6(20) (5) MeV , T ph = 184(23) MeV . 



(31) 



(32) 



where Y p h is computed using the computed value of g p7nr 
and the physical values of m p and m,. The first quoted 
uncertainty indicates the statistical uncertainty and the 
second indicates the uncertainty associated with the de- 
termination of the lattice spacing. 

In Fig. 3 we show the scattering phase shifts and 
the fit curve we obtained using the effective range for- 
mula. With the exception of the left-most and right-most 
points, the data are not more than one standard devia- 
tion away from the curve. However, the confidence level 
of the fit is Q « 8%, and it is not clear whether the ef- 
fective range formula is reliable for = 304(2) MeV. 
Assuming the effective range formula is reliable, the value 
we computed for the coupling constant is in reasonably 
good agreement with the physical value quoted by the 
PDG gP%f = 5.975(16) [15]. 

In Fig. 4 we compare our results to other recent lattice 
studies [9, 13, 21]. Our result for the coupling constant 
is in good agreement with other studies and was deter- 
mined with similar precision. The value of m p is also 
compatible with the value reported by the other studies, 
but our error bar is significantly smaller. We note that 
the results reported by Lang and collaborators [13] differ 
significantly from our results and the ones of other lattice 
studies. This is presumed to be due to the fact that the 
volume used in their study is too small. 

The phase shifts and the values we obtained for m p and 
gp-mr are also in good agreement with one-loop unitarized 
chiral perturbation theory [14]. To show the agreement, 
in the left panel of Fig. 3 we plot the phase shift expected 
from unitarized xPT to 0(p 4 ) for = 304 MeV. Note 
that this is not a fit to our data - the only input from our 
study to this expression is the pion mass. Another check 
is presented in the left panel of Fig. 4 where we plot the 
expectation for the mass of the p as a function of the 
pion mass. Our datapoint falls directly on top of the 
predicted band in the region where the unitarized xPT 



7 




(r m n f m n [MeV] 

FIG. 4: Left panel: m p versus m n from recent lattice studies. We use lattice units to remove the uncertainty associated with 
setting the scale. The band represents the prediction from unitarized x?T [14] (note that this is not a fit to our data). Right 
panel: coupling as a function of from recent lattice studies. The dot-dashed line is the PDG result [15]. 



prediction has the smallest uncertainty. This is not a 
trivial agreement since this band has no input from our 
study. 

IV. CONCLUSIONS 

We presented the first study of the p meson decay 
using asymmetrical lattices. The calculation was car- 
ried out using nHYP-smeared clover fermions with two 
mass-degenerate quark flavors at a pion mass of = 
304(2) MeV. We used three ensembles of 300 dynamical 
gauge configurations and computed the P-wave scatter- 
ing phase shifts for the isospin I — 1 channel. Fitting the 
phase shifts using the effective range formula, we found 
a coupling constant g p7r7r = 6.67(42) and resonance mass 
to p = 827(3) (5) MeV. Using the computed value for the 
coupling constant and assuming that it does not change 
significantly as the quark mass is lowered towards the 
physical point, we estimate the physical decay width to 
be Tph = 184(23) MeV, which is compatible with the 
experimental value T exp = 149.1(8) MeV. 

To compute the scattering phase shifts, we chose lat- 
tices volumes with one spatial direction elongated. By 
varying the elongation parameter, we were able to com- 
pute a set of scattering phase shifts evenly distributed 
throughout the spectral region where the p decays. While 
the moving frame formalism has been successfully used to 
study the p decay for as low as = 266 MeV on cubic 
lattices, the finer momentum control offered by asymmet- 
rical lattices will be needed to study narrower resonances. 

We found that both the Tnr(t) and p(t) interpolators 
couple to the two-pion scattering states and can be used 
to extract the ground state. However, using the varia- 
tional method allows both the ground and first excited 
state to be computed and gives a more precise determi- 
nation of the ground state. 

To compute the quark diagrams required for the two- 



pion operator, we had to employ stochastic methods. 
When employing spin-color dilution for the noise, we 
found that the gauge variance of the estimators is sub- 
stantially larger than the noise variance. It is then suffi- 
cient to use only two independent noise vectors per con- 
figuration - additional vectors will add significantly to 
the analysis cost without reducing the final error bars. 

To compute the resonances parameters, we used the 
effective range formula. While this formula described the 
overall phase shift profile well, the confidence level for our 
fit was only Q ~ 8%. We find that the effective range 
formula cannot reliably describe the phase shifts for the 
entire momentum range used in our study. However, the 
data points in the resonance region are well described by 
this formula. 

Our results arc compatible with most recent lattice 
studies and our error bars arc comparable or better. 
Moreover, our results were in very good agreement with 
the predictions of unitarized chiral perturbation theory. 

For future studies, we plan to increase the size of the 
variational basis as proposed by [10], include smeared 
quarks, and compute the resonance parameters at a pion 
mass ~ 200 MeV. The systematic effects due to lat- 
tice spacing, finite volume, and the inclusion of strange 
quark dynamics will be determined in further studies. 
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Appendix A: Numerical evaluation of Zi m (s,q ; 772, 773) 

In the following we describe the numerical evaluation 
of the generalized zeta functions defined by 



Zi m (s,q 2 ;m,V3) = X] 



Mm(n) 



(n 2 - g 2 ) 



n£2 



2\a 



for 



Re 2s < I + 3 



(Al) 



(A2) 



We use the same strategy originally proposed by Lusher 
in [3]. Its extension to asymmetrical lattices is given 
in [7]. In the following, we give the relevant formulas, 
describe the numerically procedure used in this work, and 
give a list of values for comparison. 

We begin by introducing the heat kernel 



AC(t,x) 



— Y 



mm 

(47Tt) 3 / 2 



E< 

n6Z 3 



= -3t( x - 27rft )' 



(A3) 

(A4) 



where n = (m, rj 2 n 2 , m 713) and n = (711,712/772,713/773). 
The first form is the definition of the heat kernel. 
Eq. (A4) follows from Eq. (A3) by applying the Poisson 
summation formula 



/oo 
/( 



„z27rk-x j 3 n 



(A5) 



and will be useful for small times t. The derivatives of 
the heat kernel are defined by 



K,i m (t,x.) = y;,„(-iV x )/C(t,x) 



(A6) 



where the expression 3^ m (— «V X ) is to be understood as 
the replacement of the homogeneous polynomial 



J4 m (x) = x l Yi m {n x ) , 



(A7) 



using the rule Xi — > Lastly, we define the truncated 
heat kernel and its derivatives as 

IC A (t, x) = JC(t, > ( A8 ) 



and 



/Cf m (i,x) -y /m (-iV x )/C A (t,x). (A9) 



A valid representation of the generalized zeta function, 
analytically continued to the half plane Re 2s > 1/2, can 
be written in terms of the truncated heat kernel as 

Zi m (s,q 2 ; 772,773) 



E 



^m(n) (2-n-) 3 f SioSmomm 



(n 2 -g 2 )* T(s) l(4^) 2 (s-3/2) 



<A 
1 

dtt 8 ' 1 
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(A10) 



where the cutoff A is chosen such that A 2 > Reg 2 . In 
order to compute the scattering phase shifts, we need to 
evaluate Zi m (l, q 2 ; 772, 773) numerically. In the following 
we describe the procedure used in this work. 
For the case of interest, s = 1, we have 



2i m (^,q 2 ;V2,m) 

= E 



^m(n) , , 3 f <Wm0772773 



|n|<A 
• 1 

dt 
10 



(n - q 2 ) 



Ji 2 r^ 1+ n\ <WW72?73 
e« K lm {t,0)- (47r)2t3/2 



(All) 



In order to evaluate Zi m (l, g 2 ; 77 2 , 773), the infinite sums 
appearing in K.i m (t, 0) must be truncated. In the inte- 
gration region [1, 00], we use Eq. (A3) to define K.i m (t, 0). 
In this case subsequent terms in the sum are suppressed 
by Ree-*^ 2 - 1 ? 2 ), and the sum converges the slowest for 
t = 1. In the integration region [0, 1], we use the alterna- 
tive definition of the heat kernel to define K,^ m (t,0). In 
this case one has to approximate e tq K.i m (t,0). Again, 
the sum converges the slowest at t = 1 with succes- 

(7rn) 2 2 

sive terms being suppressed by Re e t vtq . The cut- 
off is then chosen such that the integrands are approx- 
imated to a desired precisions for t — 1. In Table III 
we list some values for comparison, and in Figure 5 
we plot -2bo(l, q 2 ', 1, 2) and 2 2 o(l, q 2 ; 1, 2) in the range 
< g 2 < 4. 



-2oo(l,g 2 ,l,2) 



.220(1, g 2 , 1,2) 



0.2 

0.41 

0.62 

0.83 

1.04 



6.35553 
-5.08222 
1.16621 
10.9025 
-33.5287 



5.62348 
-2.90547 
-2.22911 
-2.85516 
-5.64013 



TABLE III: A few values for Z 00 (l, q 2 , 1, 2) and 
£20(1, g 2 , 1,2). 
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Appendix B: Determination of the lattice spacing 

To determine the lattice spacing, we use the Sommer 
scale [22]. We work in the Coulomb gauge and extract 
the heavy quark potential from the relation 



-V(r)T 



(Bl) 



which is valid for large T with Wilson line Wt defined 
by 



T-l 



W T (x,i) = Y[ C/ 4 (x,t 



i=0 



To fit the data, we use the functional form 



B / 1 

V(r) = C+ — + ar + \[- 
r \ r 



latt 



(B2) 



(B3) 



where C is part of the quarks' self-energy, a is the string 
tension, B — — \ct s , and the last term is the difference 
between the lattice and continuum potentials. The Som- 
mer parameter is then determined through 



ro/a 



1.65 + B 



(B4) 



Using a combined fit to all three ensembles and averaging 
over all space-time points, we find 

C = 0.919(4) , B = -0.430(5) , a = 0.077(1), (B5) 

and 

r Q /a = 3.984(21) . (B6) 

To fix the lattice spacing, we take r$ = 0.5 fm and obtain 
a = 0.1255(7) fm. 
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